Magnetic fluctuations and formation of large-scale inhomogeneous magnetic 

structures in a turbulent convection 
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Magnetic fluctuations generated by a tangling of tiie mean magnetic field by velocity fiuctuations 
are studied in a developed turbulent convection with large magnetic Reynolds numbers. We show 
that the energy of magnetic fluctuations depends on magnetic Reynolds number only when the 
mean magnetic field is smaller than _Beq/4Rm^'''', where Beq is the equipartition mean magnetic 
field determined by the turbulent kinetic energy and Rm is magnetic Reynolds number. Generation 
of magnetic fiuctuations in a turbulent convection with a nonzero mean magnetic field results in a 
decrease of the total turbulent pressure and may cause formation of the large-scale inhomogeneous 
magnetic structures even in an originally uniform mean magnetic field. This effect is caused by a 
negative contribution of the turbulent convection to the effective mean Lorentz force. The inhomo- 
geneous large-scale magnetic fields are formed due to the excitation of the large-scale instability. 
The energy for this instability is supplied by the small-scale turbulent convection. The discussed 
effects might be useful for understanding the origin of the solar nonuniform magnetic fields, e.g., 
sunspots. 

PACS numbers: 47.65.Md 



I. INTRODUCTION 

Magnetic fields in astrophysics are strongly nonuniform 
(see, e.g., [H, i, i, i, i, i, II i)- Large-scale magnetic 
structures are observed in the form of sunspots, solar 
coronal magnetic loops, etc. There are different mecha- 
nisms for the formation of the large-scale magnetic struc- 
tures, e.g., the magnetic buoyancy instability of stratified 
continuous magnetic field [1, [^, [13, [HI , the magnetic flux 
expulsion |12| . the topological magnetic pumping [isj . 
etc. 

Magnetic buoyancy applies in the literature for differ- 
ent situations (see [ll|). The first corresponds to the 
magnetic buoyancy instability of stratified continuous 
magnetic field (see, e.g., 0, H, IS Uli ) i and magnetic flux 
tube concept is not used there. The magnetic buoyancy 
instability of stratifled continuous magnetic field is ex- 
cited when the scale of variations of the initial magnetic 
field is less than the density stratification length. On 
the other hand, buoyancy of discrete magnetic flux tubes 
has been discussed in a number of studies in solar physics 
and astrophysics (see, e.g., [H El, El, El 113, El ) • This 
phenomenon is also related to the problem of the storage 
of magnetic fields in the overshoot layer near the bottom 
of the solar convective zone (see, e.g., [Tol. [20l. [2ll. [22j ) . 

A universal mechanism of the formation of the nonuni- 
form distribution of magnetic flux is associated with a 
magnetic flux expulsion. In particular, the expulsion of 
magnetic flux from two-dimensional flows (a single vor- 
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tex and a grid of vortices) was demonstrated in [iJI . In 
the context of solar and stellar convection, the topolog- 
ical asymmetry of stationary thermal convection plays 
very important role in the magnetic fleld dynamics. In 
particular, the topological magnetic pumping is caused 
by the topological asymmetry of the thermal convection 
The fluid rises at the centers of the convective cells 
and falls at their peripheries. The ascending fluid el- 
ements (contrary to the descending fluid elements) are 
disconnected from one another. This causes a topological 
magnetic pumping effect allowing downward transport of 
the mean horizontal magnetic field to the bottom of a cell 
but impeding its upward return 0, 

Turbulence may form inhomogeneous large-scale mag- 
netic fields due to turbulent diama gne tic a nd p aramag- 
netic effects (see, e.g., 0, HI IH M 113, H!])- Inho- 
mogeneous velocity fluctuations lead to a transport of 
mean magnetic flux from regions with high intensity of 
the velocity fluctuations. Inhomogeneous magnetic fluc- 
tuations due to the small-scale dynamo cause turbulent 
paramagnetic velocity, i.e., the magnetic flux is pushed 
into regions with high intensity of the magnetic fluctu- 
ations. Another effects are the effective drift velocities 
of the mean magnetic field caused by inhomogeneities of 
the fluid density [H, [23| and pressure [HI . In a nonlinear 
stage of the magnetic field evolution, inhomogeneities of 
the mean magnetic field contribute to the diamagnetic 
or paramagnetic drift velocities depending on the level of 
magnetic fluctuations due to the small-scale dynamo and 
level of the mean magnetic field 30]. The diamagnetic 
velocity causes a drift of the magnetic field components 
from the regions with a high intensity of the mean mag- 
netic field. 

The nonlinear drift velocities of the mean magnetic 
field in a turbulent convection have been determined in 



2 



[3l|. This study demonstrates that the nonhnear drift 
velocities are caused by the three kinds of the inhomo- 
geneities, i.e., inhomogeneous turbulence; the nonuni- 
form fluid density and the nonuniform turbulent heat 
flux. The nonlinear drift velocities of the mean magnetic 
field cause the small-scale magnetic buoyancy and mag- 
netic pumping effects in the turbulent convection. These 
phenomena are different from the large-scale magnetic 
buoyancy and magnetic pumping effects which are due to 
the effect of the mean magnetic field on the large-scale 
density stratified fluid flow. The small-scale magnetic 
buoyancy and magnetic pumping can be stronger than 
these large-scale effects when the mean magnetic field is 
smaller than the equipartition field determined by the 
turbulent kinetic energy [3ll |. The pumping of magnetic 
flux in three-dimensional compressible magnetoconvec- 
tion has been studied in direct numerical simulations in 
[s^ l by calculating the turbulent diamagnctic and para- 
magnetic velocities. 

Turbulence may affect also the Lorentz force of the 
large-scale magnetic field (see [H, HI, [H, [S^). This 
effect can also form inhomogeneous magnetic struc- 
tures. In this study a theoretical approach proposed in 
[3^ . m, [3^ for a nonconvective turbulence is further 
developed and applied to investigate the modification of 
the large-scale magnetic force by turbulent convection 
and to elucidate a mechanism of formation of inhomoge- 
neous magnetic structures. 

This paper is organized as follows. In Sect. II we dis- 
cuss the physics of the effect of turbulence on the large- 
scale Lorentz force. In Sect. Ill we formulate the gov- 
erning equations, the assumptions, the procedure of the 
derivations of the large-scale effective magnetic force in 
turbulent convection. In Sect. IV we study magnetic 
fluctuations and determine the modification of the large- 
scale effective Lorentz force by the turbulent convection. 
In Sect. V we discuss formation of the large-scale mag- 
netic inhomogeneous structures in the turbulent convec- 
tion due to excitation of the large-scale instability. Fi- 
nally, we draw conclusions in Sect. VI. In Appendix A we 
perform the derivation of the large-scale effective Lorentz 
force in the turbulent convection. 



II. TURBULENT PRESSURE AND EFFECTIVE 
MEAN MAGNETIC PRESSURE 

In this Section we discuss the physics of the effect of 
turbulence on the large-scale Lorentz force. First, let 
us examine an isotropic turbulence. The Lorentz force 
of the small-scale magnetic fluctuations can be written 
in the form J^.'™-* — ^ where the magnetic stress 



tensor cTj-™^ is given by 
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(1) 



b are the magnetic fluctuations and 5ij is the Kronecker 
tensor. Hereafter we omit the magnetic permeability 



of the fluid /x and include in the definition of 

magnetic field. The angular brackets in Eq. ([1]) de- 
note the ensemble averaging. For isotropic turbulence 
{hibj) = 5ij (b^)/3, and the magnetic stress tensor reads 
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(2) 



where = (b^) /2 is the energy density of the magnetic 
fluctuations. The magnetic pressure Pm is related to the 



magnetic stress tensor: cr. 



(m) 



-PmStj, where p„ 



Wm/i- Similarly, in an isotropic turbulence the Reynolds 



stresses {uiUj) read: (uiUj) = Sij{u )/3, and 



-pa{UiUj) 



2Wk 



(3) 



where u are the velocity fluctuations, Wk = po(u^)/2 
is the kinetic energy density of the velocity fluctuations 
and po is the fluid density. Equation ^ yields the hy- 

drodynamic pressure pv = 2Wk/'i, where u^J'' = —pv Sij. 
Therefore, the equation of state for the isotropic turbu- 
lence is given by 
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(see also [33, |3g), where p^ is the total (hydrodynamic 
plus magnetic) turbulent pressure. Similarly, the equa- 
tion of state for an anisotropic turbulence reads 
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where An = {2/3)[{u']_)/{ul) - 2] is the degree of 
anisotropy of the turbulent velocity field u = uj_ -I- it^e. 
For an isotropic three-dimensional turbulence (u^) ~ 



2(ui 



and the parameter An 

,2\ _ 



0, while for a two- 
dimensional turbulence (u^) = and the degree of 
anisotropy An 00. Here e is the vertical unit vec- 
tor perpendicular to the plane of the two-dimensional 
turbulence. 

In a two-dimensional turbulence An ^ 00 and the 
total turbulent pressure — *■ Wk- Note that the 
magnetic pressure in a two-dimensional turbulence van- 
ishes. Indeed, for isotropic magnetic fluctuations in a 
two-dimensional turbulence (bibj) — (1/2) (b^) (5^-^'', and 



therefore, cr|™^ = 

^ij ^ — ^i'j ~ 6j • 

The total energy density Wt — Wk 



(l/2)(b2)42' + = 0, where 



Wm of the ho- 
mogeneous turbulence with a mean magnetic field B is 
determined by the equation 



dWj 



It 



Wt 
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(see, e.g., |36|), where tq is the correlation time of the 
turbulent velocity field in the maximum scale Iq of turbu- 
lent motions. It is the energy source of turbulence, rj^ is 
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the turbulent magnetic diffusion and the mean magnetic 
field B is given (prescribed). The second term, Wt/tq, in 
the right hand side of Eq. ([6]) determines the dissipation 
of the turbulent energy. For a given time-independent 
source of turbulence It the solution of Eq. ([6]) is given 
by 



Wt = TO [It + 7]^iV X Bf] 



1 — exp I — 



To 
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(7) 



where Wt = WT{t ~ 0). For instance, a time- 
independent source of the turbulence exists in the Sun. 
The mean nonuniform magnetic field causes an additional 
energy source of the turbulence. In — VtC^ ^ B)^. The 
ratio In /It of these two sources of turbulence is of the 
order of 
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It 
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where is is the characteristic scale of the spatial varia- 
tions of the mean magnetic field. Since Iq <^ Lb and 



B2 < po(u 



2\(0) 



we can neglect the small magnetic 



source In of the turbulence. Thus, for t ^ tq the total 
energy density of the turbulence reaches a steady state 
Wt = const = tqIt- Therefore, the total energy den- 
sity Wt of the homogeneous turbulence is conserved (the 
dissipation is compensated by a supply of energy), i.e 



Wk + Wm — const. 



(9) 



A more rigorous derivation of Eq. ([9]) is given in Ap- 
pendix A (see Eq. (|A16|) ). Equation ([9]) implies that the 
uniform large-scale magnetic field performs no work on 
the turbulence. It can only redistribute the energy be- 
tween hydrodynamic and magnetic fluctuations. 

Combining Eqs. (H]) and ^ we can express the change 
of turbulent pressure Spr^ in terms of the change of the 
magnetic energy density SWm for an isotropic turbulence 
Sp^ = -{1/3) SWm (see (H,!!!,!!^). Therefore, the tur- 
bulent pressure is reduced when magnetic fluctuations 
are generated {SWm > 0). Similarly, for an anisotropic 
turbulence, the generation of magnetic fluctuations re- 
duces the turbulent pressure, i.e.. 



Sp^ = 
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N 



3(2 + An) 
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The total turbulent pressure is decreased also by the 
tangling of the large-scale mean magnetic field B by the 
velocity fluctuations (see, e.g., [U, SjMB' ^^'^ references 
therein). The mean magnetic field generates additional 
small-scale magnetic fluctuations due to a tangling of the 
mean magnetic field by velocity fluctuations. For a small 
energy of the mean magnetic field, B^ <C po (u^), the 
energy of magnetic fluctuations, (b^) — {h'^)^^\ caused 
by a tangling of the mean magnetic field can be written 



in the form: 



(b^) - (b 



2\(0) 



a,„(B,Rm)B2 + 0[BV(po (u'))' 



(10) 



where (b^)*^*^^ are the magnetic fluctuations with a zero 
mean magnetic field generated by a small-scale dynamo. 
Equation (fTO|) allows us to determine the variation of the 
magnetic energy SWm- Therefore, the total turbulent 
pressure reads 



Pt 



B2 



(11) 



where p^^^ is the turbulent pressure in a flow with a zero 
mean magnetic field and the coefficient qp oc am{B, Rm). 
Here we neglect the small terms ~ 0[B'*/(po (u^))^]. The 
coefficient qp is positive when magnetic fiuctuations are 
generated, and it is negative when they are damped. The 
total pressure is 



Pu 



Pk 



PB{B)^Pk 



■D2 



(12) 



where Pk is the mean fiuid pressure and Pb{B) = B^/2 
is the magnetic pressure of the mean magnetic field. Now 
we examine the part of the total pressure Ptot that de- 
pends on the mean magnetic field B, i.e., we consider 



PmiB) = Ps{B)-qpiB)^ = {l-qp)'^ 



(13) 



P + PmiB) and 



where now Ptot 
-p^^^ ■ The pressure Pm{B) is called the effective 



(see |33l 
P = Pk- 

(or combined) mean magnetic pressure. Note that both 
the hydrodynamic and magnetic fluctuations contribute 
to the combined mean magnetic pressure. However, the 
gain in the turbulent magnetic pressure Pm is not as large 
as the reduction of the turbulent hydrodynamic pressure 
Py by the mean magnetic field B. This is due to different 
coefficients multiplying by Wm and Wk in the equation 
of state (III) [see also Eq. ([5])]. Therefore, this effect is 
caused by a negative contribution of the turbulence to 
the combined mean magnetic pressure. 

We consider the case when P ^ B^/2, so that the to- 
tal pressure Ptot is always positive. Only the combined 
mean magnetic pressure PmiB) may be negative when 
qp > 1, while the pressure Pb(P) as well as the val- 
ues Pfc, Pvj Pm, Pt are always positive. When a mean 
magnetic field B is superimposed on an isotropic turbu- 
lence, the isotropy breaks down. Nevertheless Eq. (fT^]) 
remains valid, while the relationship between qp and am 
may change. 

In this Section we use the conservation law ([9]) for the 
total turbulent energy only for the elucidation of the prin- 
ciple of the effect, but we have not employed Eq. @ to de- 
velop the theory of this effect (see for details t34.,35i.i36i|). 
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In particular, the high-order closure procedure [s^, [s^ 
and the renormalization procedure |35i] have been used 
for the investigation of the nonconvective turbulence at 
large magnetic and hydrodynamic Reynolds numbers. 

In this study we investigate the modification of the 
large-scale magnetic force by turbulent convection. We 
demonstrate that the turbulent convection enhances 
modification of the effective magnetic force and causes 
a large-scale instability. This results in formation of the 
large-scale inhomogeneous magnetic structures. 



III. GOVERNING EQUATIONS AND THE 
PROCEDURE OF DERIVATION 

In order to study magnetic fluctuations and the modi- 
fication of the large-scale Lorentz force by turbulent con- 
vection we use a mean field approach in which the mag- 
netic and velocity fields, and entropy are decomposed 
into the mean and fluctuating parts, where the fluctuat- 
ing parts have zero mean values. We assume that there 
exists a separation of scales, i.e., the maximum scale of 
turbulent motions Iq is much smaller then the character- 
istic scale Lb oi the mean magnetic field variations. We 
apply here an approach which is described in (sol . [sil . [39| 
and outlined below. 

We consider a nonrotating turbulent convection with 
large Rayleigh numbers and large magnetic Reynolds 
numbers. We use the equations for fluctuations of the 
fluid velocity, u, entropy, s', and the magnetic field, b. 
The equations for velocity and entropy fluctuations are 
rewritten in the new variables v = u and s = y^po s' . 
We also use the new variable H = B/ y/po for the mean 
magnetic field, B. On the other hand, we do not use 
a new variable for magnetic fluctuations, b. Equations 
for fluctuations of fluid velocity, entropy and magnetic 
field are applied in the anelastic approximation, that is 
a combination of the Boussinesq approximation and the 
condition div (po u) = 0. The turbulent convection is re- 
garded as a small deviation from a well-mixed adiabatic 
reference state. This implies that we consider the hydro- 
static nearly isentropic basic reference state. 

Using these equations for fluctuations of fluid veloc- 
ity, entropy and magnetic field written in a Fourier space 
we derive equations for the two-point second-order cor- 
relation functions of the velocity fluctuations {vi Vj), the 
magnetic fluctuations (bibj), the entropy fluctuations 
(ss), the cross-helicity (biVj), the turbulent heat flux 
(s Vi) and (s hi). The equations for these correlation func- 
tions are given by Eqs. (|A4|) - (|A9|) in Appendix A. We 
split the tensor {hi bj) of magnetic fluctuations into non- 

helical, hij, and helical, /i-^- , parts. The helical part hlj 
depends on the magnetic helicity (see below). We also 
split all second-order correlation functions, M^-^-^\ into 
symmetric and antisymmetric parts with respect to the 
wave vector k, e.g., hij = h^j^ + hlj\ where the tensors 

^if — i^ijO^) + describes the symmetric part 



of the tensor and h\j ~ [hij(]<.) — k)]/2 determines 
the antisymmetric part of the tensor. 

The second-moment equations include the first-order 
spatial differential operators Af applied to the third- 
order moments Al'^^^^K A problem arises how to close 
the equations for the second moments, i.e., how to ex- 
press the third-order terms AfM^^^^^ through the sec- 
ond moments M^"^ (see, e.g., [H EH, \^). We wiU 
use the spectral r approximation which postulates that 
the deviations of the third-moment terms, Ai''M(^^^^(k), 
from the contributions to these terms afforded by the 
background turbulent convection, A/'M(^^^'°)(k), are ex- 
pressed through the similar deviations of the second mo- 
ments, M(^^)(k)-M(^^'0)(k): 

T{k) 

(14) 

(see, e.g., 0, H HI, SI), where T{k) is the scale- 
dependent relaxation time. In the background turbulent 
convection the mean magnetic field is zero. The r ap- 
proximation is applied for large hydrodynamic and mag- 
netic Reynolds numbers, and large Rayleigh numbers. In 
this case there is only one relaxation time r which can 
be identified with the correlation time of the turbulent 
velocity field. A justification of the r approximation for 
different situations has been performed in numerical sim- 
ulations and theoretical studies in [ii,[45|,'46, 47, 48, 4^ 
(see also review @). The r approximation is also dis- 
cussed in Sect. VI. 

We apply the spectral r approximation only for the 
nonhelical part hij of the tensor of magnetic fluctuations. 

The helical part h\j depends on the magnetic helicity, 
and it is determined by the dynamic equation which fol- 
lows from the ma gne tic helicity conservation arguments 
(see, e.g., [13, [HM IH, 111 [Hi, H^, review [| and ref- 
erences therein) . The characteristic time of evolution of 
the nonhelical part of the tensor hij is of the order of 
the turbulent time tq — lo/uo, while the relaxation time 
of the helical part of the tensor h^^'' of magnetic fluctu- 
ations is of the order of tq Rm, where Rm = IqUo/t] is 
the magnetic Reynolds number, uq is the characteristic 
turbulent velocity in the maximum scale of turbulent mo- 
tions If) and T] is the magnetic diffusivity due to electrical 
conductivity of the ffuid. 

In this study we consider an intermediate nonlinearity 
which implies that the mean magnetic field is not strong 
enough in order to affect the correlation time of turbulent 
velocity field. The theory can be expanded to the case a 
very strong mean magnetic field after taking into account 
a dependence of the correlation time of the turbulent 
velocity field on the mean magnetic field. 

We assume that the characteristic time of variation of 
the mean magnetic field B is substantially larger than 
the correlation time r(fc) for all turbulence scales. This 
allows us to get a stationary solution for the equations 
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for the second-order moments given by Eqs. (jA10|) - (|A14[) 
in Appendix A. For the integration in k-space of the sec- 
ond moments we have to specify a model for the back- 
ground turbulent convection (with a zero mean magnetic 
field, B = 0). Here we use the model of the background 
shear-free turbulent convection with a given heat flux 
(see Eqs. (jA17l) - (|A20p in Appendix). In this model ve- 
locity and magnetic fluctuations are homogeneous and 
isotropic. 

This procedure allows us to study magnetic fluctua- 
tions with a nonzero mean magnetic field and to investi- 
gate the modification of the large-scale Lorentz force by 
turbulent convection (see Sect. IV). 



IV. MAGNETIC FLUCTUATIONS AND 
LARGE-SCALE EFFECTIVE LORENTZ FORCE 

A. Magnetic fluctuations with a nonzero mean 
magnetic field 

Let us study magnetic fluctuations with a nonzero 
mean magnetic field using the approach outlined in Sect. 
III. Integration in k space in Eq. (jAlip yields an analyt- 
ical expression for the energy of magnetic fluctuations, 
(b^) [see Eq. ()A21|) in Appendix A]. The energy of mag- 
netic fluctuations versus the mean magnetic field B / B^q 
is shown in Fig. 1, where B^q is the equipartition mean 
magnetic field determined by the turbulent kinetic en- 
ergy. The asymptotic formulae for (b^) are given be- 
low. In particular, for a very weak mean magnetic field, 
B < Beq/4Rm^/^, the energy of magnetic fluctuations is 
given by 

(b2) = (b2)(") + l[(v2)(")-(b2)(0)]B2inRm 
3 



.^(v2)(0) (2-3cos2. 
5 



(15) 



where the quantities with the superscript (0) correspond 
to the background turbulent convection (with a zero 
mean magnetic field) , (v^)'^*'^ and (b2)(°^ are the velocity 
and magnetic fluctuations in the background turbulent 
convection. Here the magnetic field B is measured in the 
units of Bcq and (j) is the angle between the vertical unit 
vector e and the mean magnetic field B. The unit vector 
e is directed opposite to the gravity field. The parameter 
a* characterizing the turbulent convection is determined 
by the budget equation for the total energy, and it is 
given by 
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where is the turbulent viscosity, U is the mean fluid 
velocity and F^, = {u^ s')'"' is the vertical heat flux in the 
background turbulent convection. The energy of mag- 
netic fluctuations for a very weak mean magnetic field, 
B ^ i?cq/4Rm^/^, depends on the magnetic Reynolds 



number: (b^) oc InRm. This is an indication of that the 
spectrum of magnetic fluctuations is in the limit of 
a small yet finite mean magnetic field [H, [H, [sl, [s^l 
(see also discussion in [s^). When the mean magnetic 
field i3oq/4Rm^/^ <^ B <^ ^oq/4, the energy of magnetic 
fluctuations is given by 



(b2) = (b2)(0) + — [(v2)(")-(b2)(")]B2|ln(4B)| 

3 



+£^(v2)(0) B' (2 + 3cos2 
5 

and for B ^ i?oq/4 it is given by 



(16) 



(17) 



The normalized energy of magnetic fluctuations (b^) /B^ 
versus the mean magnetic field is shown in Fig. 2 for 
a nonconvective and convectivc turbulence. Inspection 
of Figs. 1-2 shows that turbulent convection increases 
the level of magnetic fluctuations in comparison with 
the nonconvective turbulence. It follows from Eqs. (fT5|) - 
(fTTll that in the case of Alfvenic equipartition, (u^)'*'^ = 
(b^) , a deviation of the energy of magnetic fluctuations 
from the background level is caused by the turbulent con- 
vection. 



B. The large-scale eff'ective Lorentz force 

The effective (combined) mean magnetic force which 
takes into account the effect of turbulence on magnetic 
force, can be written in the form J^f^ = S/ja°^, where 
the effective stress tensor crff reads 



~pa (uiUj) 



B,B, - - (b^) 6,, + {hh.j) 



(18) 



The last three terms in RHS of Eq. (|18p determine the 
contribution of velocity and magnetic fluctuations to the 
effective (combined) mean magnetic force. Using Eqs. 
(|A10p - (|Alip for (uiUj) and (bibj) after the integration in 
k space we arrive at the expression for the effective stress 
tensor: 



.off 



-[l-q^{B)]^S.,, 
+a,afjB) , 



[l-qs{B)]B,Bj 



(19) 



where the analytical expressions for the nonlinear coeffi- 
cients qp{B) and qs{B) are given by Eqs. (|A23p and (|A24p 
in Appendix A, the tensor a:^{B) is the anisotropic con- 
tribution caused by turbulent convection to the effective 
stress tensor (which is given by Eq. (|A22p in Appendix 
A). 
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FIG. 1: (a). The energy of magnetic fluctuations (b^) versus 
the mean magnetic field B/Baq for a nonconvective turbulence 
(a* — 0), Rm — 10^ and different values of the parameter 
e = {b2)W/(u^>'°': e = (solid line); e = 0.3 (dashed line) 
and e = 0.5 (thin dashed). 

(b). The energy of magnetic fluctuations (b"^) versus the 
horizontal (dashed line) and vertical (thin solid line) mean 
magnetic field for a convective turbulence (a, — 0.7), and for 
Rm = 10^ e = 0. 



The nonlinear coefficients qp{B) and qs{B) in Eq. ^ 
for the effective stress tensor are shown in Figs. 3a and 4a 
for different values of the magnetic Reynolds numbers. 
The nonlinear coefficients qp{B) and qs{B) increase with 
the magnetic Reynolds numbers in the range of weak 
mean magnetic fields {B < 0.1 Bgq). On the other hand, 
the turbulent convection reduces these nonlinear coeffi- 
cients in comparison with the case of a nonconvective tur- 
bulence. The asymptotic formulae for the nonlinear co- 
eSicients qp{B) and qs{B) are given below. In particular, 
for a very weak mean magnetic field, B ^ i3eq/4Rm^^'', 



FIG. 2: The normalized energy of magnetic fluctuations 
(b^)/B^ versus the mean magnetic field B/Bcq for a noncon- 
vective turbulence (a* — 0) (solid line); and for a convective 
turbulence (a* = 0.7) with the horizontal (dashed line) and 
vertical (thin solid line) mean magnetic field, where the cases 
B < Bcq (Fig. 2a) and B ~ Boq (Fig. 2b) are shown. Here 
Rm = 10^ and e = 0. 



the nonlinear coefficients qp{B) and qs{B) are given by 



qp{B) = ^(l-e) 



q,,{B) - -(1-e) 
15 



In Rm H 

45 



In Rm H 

15 



_^(ll_13cos2. 
35 



24a* 
~35" 



(20) 
(21) 



where e = (b2)(")/(u2)("). For B^jARm^^^ < B < 
-Bcq/4 these nonlinear coefficients are given by 

16 



25 



(1 - e) [5|ln(4B)| + 1 + 32B^ 



35 



(11- 13cos2 



(22) 
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FIG. 3: (a). The nonlinear coefficient qp{B) for different val- 
ues of the magnetic Reynolds numbers Rm: Rm = 10^ (thin 
soUd line); Rm = 10® (dashed-dotted line); Rm = 10^° 
(thick solid line) for a nonconvective turbulence (a* = 0), 
and at Rm = 10® (dashed line) for a convective turbulence 
(a, = 0.7). 

(b). The effective (combined) mean magnetic pressure 
Pm{B) = (1 — qp)B^/B^^ at Rm = 10® for a nonconvective 
turbulence (a, = 0) (thick solid line), and for a convective 
turbulence (o» = 0.7) for the horizontal field (dashed) and 
for vertical field (thin solid line). 




FIG. 4: (a). The nonlinear coefficient qs{B) for different val- 
ues of the magnetic Reynolds numbers Rm: Rm = 10^ (thin 
solid line); Rm = 10*' (thin dashed-dotted line); Rm = 10^° 
(thick solid) for a nonconvective turbulence a, — 0, and at 
Rm = 10® (dashed line) for a convective turbulence a* = 0.7. 
(b). The effective (combined) mean magnetic tension 
crg{B) = (1 — qs)B'^/B^q at Rm = 10® for a nonconvective 
turbulence (a* = 0) (thick solid line), and for a convective 
turbulence (a* = 0.7) for the horizontal field (dashed line) 
and for vertical field (thin solid). 



•^2 



ln(4B)| + - + 12B^ 
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and for B ^ B^q/ 4 they are given by 

1 , TTfl 

6B 



qpiB) 

qs{B) 



48B3 



80B 
37ra 



160S 



(l-3cos^ 



(23) 



(1- 5 008^0), (24) 



(25) 



The effective (combined) mean magnetic pressure 
Pm{B) = (1 — qp) B^/B'^^ is shown in Fig. 3b. Inspection 



of Fig. 3b sliows that tlio combined mean magnetic pres- 
sure Pm{B) = (1 — qp)B^/B^^ vanishes at some value 
of the mean magnetic field B = Bp ^ (0.2 — 0.3) Bcq. 
This causes the following effect. Let us consider an iso- 
lated tube of magnetic field lines. When B = Bp, the 
combined mean magnetic pressure Pm{B) = 0, the fluid 
pressure and fluid density inside and outside the isolated 
tube are the same, and therefore, this isolated tube is in 
equilibrium. When B > Bp, the combined mean mag- 
netic pressure P„i{B) > 0, the fluid pressure and fluid 
density inside the isolated tube are smaller than the fluid 
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pressure and fluid density outside the isolated tube. This 
results in upwards floating of the isolated tube. On the 
other hand, when B < Bp, the combined mean magnetic 
pressure P,n{B) < 0, the fluid pressure and fluid density 
inside the isolated tube are larger than the fluid pressure 
and fluid density outside the isolated tube. Therefore, 
this isolated magnetic tube flows down. 

The effective (combined) mean magnetic tension 
a^{B) = (f - qs)B^/B^^ is shown in Fig. 4b. The com- 
bined mean magnetic tension (B) vanishes at some 
value of the mean magnetic field B — Bg ^ 0.2 Bf.^ (see 
Fig. 4b). When B > Bs, the combined mean magnetic 
tension cr^ (B) > 0, and Alfvenic and magneto-sound 
waves can propagate in the isolated tubes. On the other 
hand, when B < Bs, the combined mean magnetic ten- 
sion CT^ (B) < 0, and Alfvenic and magneto-sound waves 
cannot propagate in the isolated tubes. 

The anisotropic contributions a{j{B) to the effective 
stress tensor determine the anisotropic mean magnetic 
tension due to the turbulent convection. The ten- 
sor cr^{B) is characterized by the function a^{B) = 
af'j{B)eij = QeB^/B^^. The nonlinear coefficient qe{B) 
and the anisotropic mean magnetic tension, (B) are 
shown in Fig. 5a and 5b. In the next Section we show that 
the anisotropic mean magnetic tension (B) caused by 
the turbulent convection, strongly affects the dynamics 
of the horizontal mean magnetic field. 

In this Section we demonstrate that turbulent convec- 
tion strongly modifies the large-scale magnetic force. Let 
us discuss a possibility for a study of the effect of tur- 
bulence on the effective (combined) mean Lorentz force 
in the direct numerical simulations. Consider the mean 
magnetic field which is directed along the a;-axis, i.e., 
B = Be^. Let us introduce the functions (Tx{B) and 
ay{B): 

^.{B) = -\{h^) + {hi) - p^{ul) , (26) 
^y{B) = -\{h') + {hl)-p,{ul), (27) 

which allow us to determine the coefficients qp{B) and 
qs{B) in the effective stress tensor: 

IpiB) = ^WyiB) - <7y{B = 0)] , (28) 

qsiB) = ^qp{B) - -^[a,(B) - <J,{B ^ 0)] . (29) 

Therefore, Eqs. allows to determine the effec- 

tive (combined) mean magnetic force in the direct nu- 
merical simulations. 

V. THE LARGE-SCALE INSTABILITY 

The modification of the mean magnetic force by the 
turbulent convection causes a large-scale instability. In 
this study we investigate the large-scale instability of 




eq 



FIG. 5: The nonlinear coefficient qe{B) versus the mean mag- 
netic field - (a) and the anisotropic mean magnetic tension, 
{B) = (jf^Bij = qeB^ /B^q - (b) for a convective turbulence 
(a* = 1) with the horizontal mean magnetic field (solid line) 
and vertical mean magnetic field (dashed line). 



continuous magnetic field in small-scale turbulent con- 
vection and we do not consider buoyancy of the discrete 
magnetic flux tubes. In order to study the large-scale 
instability in a small-scale turbulent convection we use 
the equation of motion (with the effective magnetic force 
!F°^ = Vjcr^j^ determined in Section IV), the induction 
equation and the equation for the evolution of the mean 
entropy [see Eqs. I|A25|) - (|A27P in Appendix A]. We es- 
timate the growth rate 7 of this instability and the fre- 
quencies uj of generated modes neglecting turbulent dis- 
sipative processes for simplicity's sake. We also neglect 
very small Brunt- Vaisala frequency based on the gradi- 
ent of the mean entropy. We seek for the solution of these 
equations in the form cx exp{-^t + icut — iK • R). 
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A. Horizontal mean magnetic field 



First, we study the large-scale instability of a horiz 
tal mean magnetic field that is perpendicular to the gr 
ity field. Let the z axis of a Cartesian coordinate syst 
be directed opposite to the gravitational field, and let 
X axis lie along the mean magnetic field B. Consider 
case = that corresponds to the interchange mo 
The dispersive relation for the instability reads 



Ky Ca 

KLr. 



(Q + iD) 



( 



where a — j -\- iuj, Lp is the density stratification leng 
Lg is the characteristic scale of the mean magnetic fi 
variations, Ca = B / ^/p^ is the Alfven speed and 



L 



+ (1-4 



2/<(2/) 1-2 



p_ 

Lb 



<^a{v) = 'le{y)y, cr'^iy) = da^{y)/dy, B is measured in 
the units of the equipartition field Bg^ = ^/pqUq and 



When Q > the growth rate of perturbations with the 
frequency 



to 



is given by 



Ca 



D 



D-^ + Q 



-1/2 



7 = 



1/2 



(31) 



(32) 



When Q < the growth rate of perturbations with the 
frequency 



UJ 



-sgn{D) 



Ca 



IQI)'^',(33) 



is given by 



7 ^ 



D^ + \Q\ 



-1/2 



(34) 



Therefore, in small-scale turbulent convection this large- 
scale instability causes excitation of oscillatory modes 
with growing amplitude. In a nonconvective turbulence 
(a* = 0) this large-scale instability is aperiodic. 

The growth rate 7 of the large-scale instability and fre- 
quency LD of the generated modes for the horizontal mean 
magnetic field versus B/B^q for different values and 
Lpl Lb for a nonconvective and convective turbulence are 
shown in Figs. 6 and 7. Here 7 and to are measured in 




0.2 0.4 0.6 s/s, 
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FIG. 6: The growth rate (a) of the large-scale instability and 
frequency of the generated modes (b) for the horizontal mean 
magnetic field with Lp/Lb ~ 0.2 versus B/Bcq for different 
values K^: Lp = 2, (dashed-dotted line) and Kz Lp = 5 
(solid line) for turbulent convection with a* = 1 (thick curves) 
and a, = 0.3 (thin curves), and for a nonconvective turbulence 
a* = (dashed line). Here e = and Rm — lO'' . 



the units of t;7^, where = {Lp/uo){K/Ky). In the tur- 
bulent convection there are two ranges for the large-scale 
instability of the horizontal mean magnetic field (when 
Q > and Q < 0), while in a nonconvective turbulence 
(a* — 0) there is only one range for the instability. The 
first range for the instability is related to the negative 
contribution of turbulence to the effective magnetic pres- 
sure for the case of Lp < Lb, while the second range is 
mainly caused by the anisotropic contribution (B) due 
to turbulent convection. 

In the absence of turbulence (small Reynolds num- 
bers) or turbulent convection (small Rayleigh numbers) 
the coefficients Q — Lp/LB — 1, D — 0, and the crite- 
rion for the large-scale instability, Lp > Lb, coincides 
with that of the Parker's magnetic buoyancy instabil- 
ity (0, S this case uj — 0, i.e., the oscillatory 
modes with growing amplitude are not excited. On the 
other hand, in a developed turbulent convection the ef- 
fective (combined) magnetic pressure becomes negative 
and the Parker's magnetic buoyancy instability cannot 
be excited. However, the instability due to the modifica- 
tion of the mean magnetic force by small-scale turbulent 
convection can be excited even when Lp < Lb, i.e., even 
in uniform mean magnetic field. 

The instability mechanism due to the modification of 
the mean magnetic force consists in the following. An iso- 
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FIG. 7: The growth rate (a) of the large-scale instability and 
frequency of the generated modes (b) for the horizontal mean 
magnetic field with Lp/Ls — 10 versus B/Beq for different 
values Kz'. K^Lp = 3 (dashed-dotted line) and K^Lp — 5 
(solid line) for turbulent convection with a, = 1 (thick curves) 
and a, = 0.3 (thin curves), and for a nonconvective turbulence 
a* = (dashed line). Here e = and Rm = lO''. 



lated tube of magnetic field lines moving upwards, turns 
out to be lighter than the surrounding plasma. This is 
due to the fact that the decrease of the magnetic field 
in the isolated tube caused by its expansion, is accom- 
panied by an increase of the effective magnetic pressure 
inside the tube. Since the effective (combined) magnetic 
pressure is negative, this leads a decrease of the density 
inside the tube. The arising buoyant force results in the 
upwards floating of the isolated tube, i.e. it causes the 
excitation of the large-scale instabihty (see also discus- 
sion in |i36l] ) . 



B. Vertical uniform mean magnetic field 

Consider vertical uniform mean magnetic field, i.e., the 
magnetic field is directed along z axis. The growth rate 
of perturbations is given by 



7 = CaK, 



1/2 



(35) 



where Kj_ is the component of the wave vector that is 
perpendicular to z axis and K = ^ + . It follows 
from Eq. that the large-scale instability of the verti- 
cal uniform mean magnetic field is caused by the modifi- 



cq) 



cation of the mean magnetic tension by small-scale tur- 
bulent convection. When (?s > 1 (i.e., when B < 0.2B, 
the instability occurs for an arbitrary value K± , while for 
Qs < 1 the necessary condition for the instability reads 
K± > K where 



1 - <ls{y) 



22/k^(j/)l 



y=B^ 



and we take into account that q'g{y) < 0. The growth rate 
of the instability is reduced by the turbulent dissipation 
oc 1/^ . The maximum growth rate 7inax at a fixed 
value of the wave number K (i.e., at a fixed value of 
the turbulent dissipation) is attained at = = 
K [(1 - x)/2]i/2, and it is given by 



7n 



KCa{1-x) 



yWs{y)\Y'^ 

y=B^ 



(36) 



where x < 1. The maximum growth rate 7max of the 
instability for the vertical uniform mean magnetic field 
versus B/ Bcq is plotted in Fig. 8. Here 7niax is measured 
in the units of uq K. The value of 7niax is larger for a non- 
convective turbulence, but the range for the instability is 
wider for the turbulent convection (see Fig. 8). 



VI. DISCUSSION 

In the present study we investigate magnetic fluctua- 
tions generated by a tangling of the mean magnetic field 
in a developed turbulent convection. When the mean 
magnetic field B <^ Beq/4Rm^^'*, the energy of magnetic 
fluctuations depends on magnetic Reynolds number. We 
study the modification of the large-scale magnetic force 
by turbulent convection. We show that the generation of 
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magnetic fluctuations in a turbulent convection results in 
a decrease of the total turbulent pressure and may cause 
formation of the large-scale magnetic structures even in 
an originally uniform mean magnetic field. This phe- 
nomenon is due to a negative contribution of the turbu- 
lent convection to the effective mean magnetic force. 

The large-scale instability causes the formation of in- 
homogeneous magnetic structures. The energy for these 
processes is supplied by the small-scale turbulent con- 
vection, and this effect can develop even in an initially 
uniform magnetic field. In contrast, the Parker's mag- 
netic buoyancy instability is excited when the density 
stratification scale is larger than the characteristic scale 
of the mean magnetic field variations (see The 
free energy in the Parker's magnetic buoyancy instability 
is drawn from the gravitational field. The characteristic 
time of the large-scale instability is of the order of the 
Alfven time based on the large-scale magnetic field. 

We study an initial stage of formation of the large-scale 
magnetic structures for horizontal and vertical mean 
magnetic fields relative to the vertical direction of the 
gravity field. In the turbulent convection there are two 
ranges for the large-scale instability of the horizontal 
mean magnetic field. The first range for the instability is 
related to the negative contribution of turbulence to the 
effective magnetic pressure for the case of Lp < Lb, while 
the second range for the instability is mainly caused by 
the anisotropic contribution of the turbulent convection 
to the effective magnetic force. The large-scale instability 
of the vertical uniform mean magnetic field is caused by 
the modification of the mean magnetic tension by small- 
scale turbulent convection. The discussed effects in the 
present study might be useful for the understanding of 
the origin of the sunspot formation. 

Since in the present study we neglect very small Brunt- 
Vaisala frequency based on the gradient of the mean en- 
tropy, we do not investigate the large-scale dynamics of 
the mean entropy. This problem was addressed in (soj 
whereby the modification of the mean magnetic force by 
the turbulent convection was not taken into account. 

In order to study magnetic fiuctuations and the mod- 
ification of the large-scale Lorentz force by turbulent 
convection we apply the spectral r approximation (see 
Sect. III). The r approach is an universal tool in tur- 
bulent transport that allows to obtain closed results and 
compare them with the results of laboratory experiments, 
observations and numerical simulations. The r approx- 
imation reproduces many well-known phenomena found 
by other methods in turbulent transport of particles and 
magnetic fields, in turbulent convection and stably strat- 
ified turbulent fiows (see below). 

In turbulent transport, the r approximation yields cor- 
rect formulae for turbulent diffusion, turbulent thermal 
diffusion and turbulent barodiffusion (see, e.g., [60l.l6l|). 
The phenomenon of turbulent thermal diffusion (a non- 
diffusive streaming of particles in the direction of the 
mean heat fiux), has been predicted using the stochastic 
calculus (the path integral approach) and the r approx- 



imation. This phenomenon has been already detected 
in laboratory experiments in oscillating grids turbulence 
[6^ and in a multi-fan turbulence generator in sta- 
bly and unstably stratified fluid fiows. The experimental 
results obtained in [63, l62] are in a good agreement with 
the theoretical studies performed by means of different 
approaches (see [H^ilil]). 

The T approximation reproduces the well-known 
/c~^/'^-spectrum of anisotropic velocity fluctuations in a 
sheared turbulence (see j65i|). This spectrum was found 
previously in analytical, numerical, laboratory studies 
and was observed in the atmospheric turbulence (see, 
e-g-, [11]). In the turbulent boundary layer problems, the 
T-approximation yields correct expressions for turbulent 
viscosity, turbulent thermal conductivity and the classi- 
cal heat fiux. This approach also describes the counter 
wind heat flux and the Deardorff's heat fiux in convec- 
tive boundary layers (see [65|). These phenomena have 
been studied previously using different approaches (see, 
e.g., |ill,E3,[^). 

The theory of turbulent convection [6^1 based on the r- 
approximation explains the recently discovered hystere- 
sis phenomenon in laboratory turbulent convection [68| . 
The results obtained using the r-approximation allow 
also to explain the most pronounced features of typical 
semi-organized coherent structures observed in the at- 
mospheric convective boundary layers (" cloud cells" and 
"cloud streets") The theory [H based on the r- 

approximation predicts realistic values of the following 
parameters: the aspect ratios of structures, the ratios 
of the minimum size of the semi-organized structures to 
the maximum scale of turbulent motions and the char- 
acteristic lifetime of the semi-organized structures. The 
theory also predicts excitation of convective-shear 
waves propagating perpendicular to the convective rolls 
("cloud streets"). This waves have been observed in 
the atn iosp heric convective boundary layers with cloud 
streets |69| . 

A theory [t^ for stably stratified atmospheric turbu- 
lent flows based on the r-approximation and the budget 
equations for the key second moments, turbulent kinetic 
and potential energies and vertical turbulent fluxes of 
momentum and buoyancy, is in a good agrement with 
data from atmospheric and laboratory experiments, di- 
rect numerical simulations and large-eddy simulations 
(see detailed comparison in Sect. 5 of [70|). 

The detailed verification of the r approximation in the 
direct numerical simulations of turbulent transport of 
passive scalar has been recently performed in [46'|. In 
particular, the results on turbulent transport of passive 
scalar obtained using direct numerical simulations of ho- 
mogeneous isotropic turbulence have been compared with 
that obtained using a closure model based on the r ap- 
proximation. The numerical and analytical results are in 
a good agreement. 

In magnetohydrodynamics, the r approximation repro- 
duces many well-known phenomena found by different 
methods, e.g., the r approximation yields correct formu- 
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lae for the a-effect [1, [H, [7l|, , the turbulent diamag- 
netic and paramagnetic velocities [3, EE, [H, [13, [Hi , the 



SlsflQ^the r2xJ 

the shear-current effect 



turbulent magnetic difFusion [ 
effect and the /t-effect 0, [li 

Generation of the large-scale magnetic field in a non- 
helical turbulence with an imposed mean velocity shear 
has been recently investigated in ^75i] using direct numer- 
ical simulations. The results of these numerical simula- 
tions are in a good agreement with the theoretical pre- 
dictions based on the r approximation (see [sol [39l U^) 
and with the numerical solutions of the nonlinear dynamo 
equations performed in (tgI . Ftt} (see detailed comparison 
in Hi). 

The validity of the r approximation has been tested in 
the context of dynamo theory, in direct numerical sim- 
ulations in [13]. The alpha effect in mean field dynamo 
theory becomes proportional to a relaxation time scale 
multiplied by the difference between kinetic and current 
helicities. It is shown in [i^l that the value of the relax- 
ation time is positive and, in units of the turnover time 
at the forcing wavenumber, it is of the order of unity. Ki- 
netic and current helicities are shown in (47j | to be domi- 
nated by large scale properties of the flow. Recent stud- 
ies in [4^ of the nonlinear alpha effect showed that in 
the limit of small magnetic and hydrodynamic Reynolds 
numbers, both the second order correlation approxima- 
tion (or first-order smoothing approximation) and the r 
approximation give identical results. This is also sup- 
ported by simulations ^] of isotropically forced helical 
turbulence whereby the contributions to kinetic and mag- 
netic alpha effects are computed. The study performed in 
[4^ provides an extra piece of evidence that the r approx- 
imation is a useable formalism for describing simulation 
data and for predicting the behavior in situations that 
are not yet accessible to direct numerical simulations. 
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APPENDIX A: VELOCITY AND MAGNETIC 
FLUCTUATIONS IN TURBULENT 
CONVECTION 

In order to study the velocity and magnetic fluctua- 
tions with a nonzero mean magnetic field and to derive 
the effective stress tensor in the turbulent convection, we 
use a mean field approach in which the magnetic and ve- 
locity fields, and entropy are decomposed into the mean 
and fluctuating parts, where the fluctuating parts have 



zero mean values. The equations for fluctuations of the 
fluid velocity, entropy and the magnetic fleld are given 
by 



1 av(x,i) 



Pa 



g ^ 1 

s + 



Po 

An 



(b- V)H 



-(H.V)b+ ^[2e(b.H) 



-(b.e)H] 



(Al) 



9b(x, t) 
dt 



9s(x, t) 
dt 



(H- V)v- (v V)H- 
-H(v•e)]+b^, 

-^(ve) + .^, 
9 



[v(H . e) 



(A2) 
(A3) 



where we used new variables (v, s, H) for fluctuating 
flelds V = and s = ^fp^s' , and also for the mean 

field H ~ B/^/po. Here B is the mean magnetic field, 
Po is the fiuid density, e is the vertical unit vector di- 
rected opposite to the gravity field, fi^ = — g • VS' is the 
Brunt- Vaisala frequency, S is the mean entropy, g is the 
acceleration of gravity, u, b and s' are fluctuations of 
velocity, magnetic field and entropy (we have not used 
new variables for magnetic fluctuations), v^, and 
are the nonlinear terms which include the molecular 
viscous and diffusion terms, p = p' -\- ^/po (H • b) are the 
fluctuations of total pressure, p' are the fluctuations of 
fluid pressure. 

Equations (IAip - (IA3P for fluctuations of fluid velocity, 
entropy and magnetic field are written in the anelastic 
approximation, which is a combination of the Boussi- 
nesq approximation and the condition div (po u) = 0. 
The equation, divu = Ap(u • e), in the new variables 
reads: divv — (Ap/2)(v • e), where Vpo/po = —Ape. 
The quantities with the subscript "0" correspond to the 
hydrostatic nearly isentropic basic reference state, i.e., 
VPo = Po g and g • [(7Po)^^ VPo - Po^^Po] « 0, where 
7 is the specific heats ratio and Po is the fluid pressure in 
the basic reference state. The turbulent convection is re- 
garded as a small deviation from a well-mixed adiabatic 
reference state. 

Using Eqs. l|Al |) -(|A3 p and performing the procedure 
described in Section III we derive equations for the 
two-point second-order correlation functions of the ve- 
locity fluctuations fij = (viVj), the magnetic fluctua- 
tions hij — {bibj), the entropy fluctuations — (ss), 
the cross-helicity gij ~ (biVj), the turbulent heat flux 
Fi — (svi) and Gi — {sbi). The equations for these 
correlation functions are given by 



dt 
dhjj (k) 
dt 

dt 



i(k.H)$,j +A^/„- , 
-i(k-H)$„- +Nh^j , 
i(k.H)[/,,(k)-;^,,(k) 



(A4) 
(A5) 
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dt 



aGz(k) 

dt 
5e(k) 

dt 



+genPjn{k)G,{-k)+Mg,j , (A6) 
-i(k.H)G,(k) +5e„P„,(A;)e(k) +A^F, , 

(A7) 

-i(k-H)F,(k)+A^G, , (A8) 



derivatives. Equations (jA10[) and (|A11[) yield 



1 + 2V' 



(°^)(k) 

r5e„P„(fc)^^|'^(k) . (A15) 



Therefore, when the mean heat flux F.j;^'^^ in the back- 
ground turbulence is zero (i.e., for the nonconvective tur- 
bulence), we obtain 



-^F,(k)+AAe, 

.9 



(A9) 



(see for details [311), where <&ij(k) = 3ij(k) — 5ji(— k), 
■^/zj - .ge„[/^„(fc)Fj(k) + Pj„(fc)F,(-k)] + A^/,,, and 
AA/y, TVA/iij, A/'ffy, AfFi, AfGi and A/'O are the third- 
order moment terms appearing due to the nonlinear 
terms. The terms ^ Fi in the tensor JVfij can be consid- 
ered as a stirring force for the turbulent convection. Note 
that a stirring force in the Navier-Stokes turbulence is an 
external parameter. 

We split the tensor of magnetic fluctuations into non- 
helical, hij, and helical, hlj , parts. The helical part 

depends on the magnetic helicity, and it is de- 
termined by the dynamic equation which follows from 
the magnetic helicity conservation arguments. We also 
split all second-order correlation functions into symmet- 
ric and antisymmetric parts with respect to the wave 



/^;^(k) 



-V . v^v .V . '(k) . (A16) 

This is in agreement with the fact that a uniform mean 
magnetic field performs no work on the turbulence (with- 
out mean heat flux). It can only redistribute the energy 
between hydrodynamic fluctuations and magnetic fluc- 
tuation. A change of the total energy of fluctuations is 
caused by a nonuniform mean magnetic field. For the in- 
tegration in k-space of these second moments we have to 
specify a model for the background turbulent convection 
(with zero mean magnetic field, B = 0). Here we use the 
following model of the background turbulent convection 
[denoted with the superscript (0)] : 

= po(u2)(")W^(fc)P,,(k) , 



As) r{a) 
Jij Jij : 



where the tensors 



vector k, e.g., /j, 

•l^if ~ [/«i(k) + /ii(^k)]/2 describes the symmetric part 

„(a) 



eW(k) 



= (b2)(")W^(fc)P,,(k), 
= 3po(u»sV°) VK(fc)e,Py(k) 



= 2/90 (s 



'2\(0) 



W(k) 



of the tensor and = [/,j(k) - /jj(-k)]/2 determines where Py (k) = % - hkj/k'^, W{k) 



the antisymmetric part of the tensor. We use the spectral 
T approximation (see Eq. in Sect. III). We assume 
also that the characteristic time of variation of the mean 
magnetic field B is substantially larger than the corre- 
lation time T(fc) for all turbulence scales. This allows 
us to get a stationary solution for the equations for the 
second-order moments: 



(A17) 
(A18) 
(A19) 
(A20) 

E{k)/8Trk^, 



T{k) = 2Tof(fc), E{k) = -df{k)/dk, f{k) = (fc/fco)i-9, 
1 < g < 3 is the exponent of the kinetic energy spectrum 
(e.g., (7 = 5/3 for Kolmogorov spectrum), fco — 1/Iq and 



4;'(k) 



1 



[(l + ^)/f)(k)-|-^/.(f (k) 



To = h/uQ. Note also that ^.^(k) = and G^'ilC) = 0. 

This procedure allows to study magnetic fluctuations 
with a nonzero mean magnetic field and to derive the 
effective stress tensor in the turbulent convection (see 
Section III). In particular, integration in k space in 
Eq. (|A11[) yields the energy of magnetic fluctuations 



2^Tge^P^n{k)F]''m , (AlO) 
[^/f'(k) + (l + ^)/.if)(k) 

+7AT5e„P,„(fc)Pf '(k)] , (All) 

- ^[/f^'(k)-/.r(k) 



(0) 



.2\(o)^ 



6-34"^(4P) 



4°)(4P)] 



r2\(0) 



-(l + 3cos2 (j))-^{A2} 



(A21) 



-Tge„P,„(fc)pj^^(k)] 
f^^(k) 



'''^ - l-,,/2' 
G("'(k) « -zr(k.H)i^(^)(k), 



(A12) 
(A13) 
(A14) 



(see for details [3l|), where ■!/'(k) = 2(Tk-H)^ and we 
neglected terms 0(^1). In Eqs. (|A10p - (|A14p we ne- 
glected also the large-scale spatial derivatives. The cor- 
relation functions f'>'^\ h^ij\9\f ^/"^ and c'f'' vanish 
because they are proportional to the first-order spatial 



where ^{X} = X'^^\2B) - X^^\AB), 4> is the an- 
gle between the vertical unit vector e and the mean 
magnetic field B, the functions An\y) are given by 

Eqs. (|M5|-(|X36|) and the functions A^n\v) are 

given by Eq. (IA34p . The magnetic stress tensor is given 
by Eqs. p^ . where the anisotropic contribution af^ to 
the magnetic stress tensor is determined by 

1 



ey*{2G|^^ + A^^^ + 4^'} + cos (j){ei(}j 



+e,A)^{2Gf)-4^^} 



(A22) 
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and the nonlinear coefficients qp{B) and qs{B) are given 

by 



12^2 



(0), 



(1) 



+ cos^0*{6Cf^ +4'^}] 



9.(5) = - 



(l-e)4°)(4i?) + 6a4*{C^^n 



12S2 
cos2 0*{cf^}] 



(A24) 



The asymptotic formulas for these coefficients are given 
by Eqs. ®-(I2Sl). 

In order to study the large-scale instability we use the 
equation of motion, the induction equation and the equa- 
tion for the mean entropy: 



DU, 



Dt 



= -V, 





1 

H 




Po . 



■g2 

{^-%{B))—Kpe, 



B • V)[(l - qs{B))B,] + V,[2pQV^{B) {dU),j] 

-g5-i/^(B)Ape,divU, (A25) 

(A26) 
(A27) 



— = Vx[UxB-7?,(i?) (VxB)] , 



where 77^ [B) is the turbulent magnetic diffusion, D / Dt = 
d/dt+(U-V), and F^^) = -k|J^(B)VS' is the turbulent 

(T) 

heat flux, is the tensor for the nonlinear turbulent 
thermal diffusivity, and 

Aot = Pfe + (l-gp(S)) — -1/^(5) podivU, 

Pfc is the mean fluid pressure, 2 {dU)ij — C/j Vj Ui 
and (_B) is the turbulent viscosity. The turbulent vis- 
cosity j/y (B) in Eq. (|A25p is given by 



2A[''^ + (1 + e)A'2' - 2(1 - e)H{Ai} 
(1 - 29e)cf ^ - 4(7 - 8e)H{Ci} 
(1 - 3e)G{C7i} - 2(1 - e)Q{Ci}' 



(A28) 

The explicit form of the functions H{X}, G{X} and 
Q{X} is given in [s^l- The asymptotic formula for v^{B) 
for a weak mean magnetic field, B ^ i3cq/4, is given 
by Vt{B) = v^{\ + 2e), and for B ^ i3eq/4 it is given 
by i^r(-S) = [vt/^B) (1 -I- e). The turbulent heat flux 



F^'^^ (B) = — KiJ'' (B)VjS' and the tensor for the nonlinear 
turbulent thermal diffusivity: 



«!P(B) 



(A23) where kI"^^ — /3y = BiBj/B'^. The asymptotic 



[(24°^(2B)-f 4°^(2B))(5, 



4°) (25) A,], 



(A29) 



formula for (B) for B < Beq/2Rm^/'' reads 



-!P(B) 
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■ [2(10 - f3^ In Rm) (5^ + /J^ In Rm 



(A30) 



where /3 = VSB/B^q. When Beq/2Rm^/^ < B < 
Bcq/2, the hmction «:-P(B) is 



-!P(B) 



-[(5-2/32|ln/3|)5,,-H/32|ln/3|A,] , 



and when B ^ Bcq/2, it is 



4'^(B) = .P^(%+/3.) 



(A31) 



(A32) 



In order to integrate in Eqs. (jA10[) - (jA14[) over the angles 
in k-space we used the following identity: 



ijmn 



5- d9 dip ^ Ai5^j + A2P^J , 

1 + a cos"^ 

hjmnSme ^^^^ 



1 + a cos^ 

mn 

jnf^im ~^ ^rnnPij'} i 



mn 
'in 



Hi 



(a) 



kijrnn Sill 

(1 -f acos^ 0)2 

a 

9a 



(i0 dip 



^ijmn (^) ^ q -^ijmn ('^) 



Gijmn (^) 



^iimn sin 



(1 acos2 6')3 
a 9 
2ftz 



Hijmn ('^) (-) c) Bijrnn (^^) 



where a = P'^/T{k), and 



^1 

C2 
C3 



arctan(V^) 
(a + 1) ^ 1 

^^^3^arctan(75)_^- 



27r 
a 

27r 
a 

TT 
2^^2 

A2 - 7Ai -1- 35Ci , 
Ai - 5C\ . 



(a+1) 



_^^2 3'i'ctan(Y^) 5a _^ 
\/a 3 
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(for details, see [lOl)- The functions Ail^\p) are given 
by 



where /3 = VSB/B^q. For B < Bcq/4Rm^/'' these func- 
tions are given by 



A(^HP) = ^ r^^'^^pdX, (A33) 



, (A34) 



Af\(3) - 2- -/32lnRm 
5 



2 , 



In Rm H 
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and similarly for ci™^(/3), where = P'^{k/ko)'^^'^ = a. 
The explicit form of the functions and Cii™^(/3) 

for TO = 1; 2 are given in [sO], and the functions A^^\f3) 
and Aj"'' are given by 



arctan/3, 6 o 

2 + 2 3 + 5/3^ _ /32 inRm 

Vl + /32^/R^^/ 

arctan/3,„ 9 ^9 , -r, 

2 ^ (9 + 5/32) _^ _ _ ^2 1^ 



1 + /3VRm/ 



(A36) 



For _Boq/4Rm^^'* <^ B <^ i?cq/4 these functions are given 

by 



21n/3-l| + V 
15 7 



5 



41n/5---3/32 
15 



(A35) and for B > Beq/4 they are given by 



TT 6 
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